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1. Introduction 

The effect of optical bistability (or multistability) is a basis of numerous applications 
such as optical switching, differential amplification, unidirectional transmission, power 
limiting, pulse shaping, and optical digital data processing [T]. A common way to 
obtain optically bistable devices is by using a resonant cavity containing a Kerr-type 
or saturable absorption material whose complex refraction index depends on the input 
light intensity. In the first case, the dispersion bistability appears as a result of the 
nonlinear increasing or decreasing of the optical path length as the field intensity 
inside the resonator rises. In the second one, the dependence of power absorption 
of the material inside the cavity on the input field intensity leads to the saturable 
absorption bistability. In the present paper we mainly focus our attention on the 
latter phenomenon. 

The mechanism of saturable absorption is related to the possibility of a material 
absorbing more than one photon before relaxing to the ground state. Also, two 
photons and multiphotons absorptions are known under the sufficiently high light 
intensity. In addition, population redistribution induced by intense laser fields leads 
to interaction of stimulated emission and absorption in complex molecular systems, 
and the generation of free carriers in solids [3J [3] . These phenomena are manifested 
in the both reduced (saturable) and increased (reverse saturable) optical absorption. 

The basic principle of obtaining optical bistable switching consists in using a 
saturable absorber to detune a resonant cavity by variation of the input light intensity 
[I], [4]. At low intensity, the saturable absorber is characterized by high attenuation 
which detunes the cavity and, thereby, causes the cavity to refiect and absorb 
substantially all the incident wave power. At an intensity above a threshold level of 
the system, the attenuation decreases abruptly, and the cavity transmits substantially 
all the incident wave energy. As light intensity is decreased, the cavity continues to 
transmit light until a lower threshold level is reached, at which point the medium 
again becomes a high absorber and the cavity once again is detuned. The described 
process is characterized by dependence between the input and output intensities in 
the form of a hysteresis loop. 

The threshold levels and the total change in attenuation depend upon a kind of the 
material. The key parameters for a saturable absorption material are its wavelength 
range (where it absorbs), its dynamic response (how fast it recovers), and its saturation 
intensity and fiuence Semiconductor materials, combined with proper epitaxial 
growth and correct optical design of the structure, can achieve a broad range of desired 
properties for nearly ideal saturable absorber structures for all solid-state lasers. Also 
the glasses doped by rare-earth elements or semiconductor quantum dots (QD's) have 
also been proposed as saturable absorbers to obtain optical switching in the IR range 

mm- 

Typically a saturable absorber is integrated inside a Fabry-Perot structure or 
an optical ring cavity [HI [9]. In the further realizations the intracavity saturation 
absorber is integrated in a more general mirror structure that allows for both 
saturable absorption and negative dispersion control, which is now generally referred 
to as semiconductor saturable absorber mirrors (SESAM's) [10]. In a general 
sense the design problem of SESAM's is reduced to the analysis of multilayered 
interference filters for a desired nonlinear reflectivity response. But in any case, some 
resonant cavity is used to provide the field confinement and the nonlinear response 
enhancement. 
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A perspective way to reduce the cavity size and the input field intensity required 
for optical bistable switching lies in using plasmonic and metamaterial technologies 
. It is known that the planar metamaterials can create an environment equivalent 
to a resonant cavity which allows us to achieve the required field confinement within 
a system to provide the nonlinear effect enhancement. To date there are a set of 
pubhcations [T2j [131 [TU [HI [HI [III [18] where the properties of nonlinear metamaterials 
formed by integrating nonlinear components or materials into the metamaterial are 
studied both theoretically and experimentally. Typically such structures are composed 
of metallic inclusions in the form of symmetrical split-ring resonators which are 
resonant because of an internal capacitance and inductance within each element. 
These elements are made nonlinear and tunable via the insertion of diodes with a 
voltage-controlled characteristic in the capacitive gaps of the metamaterial elements. 
To the best of our knowledge, in this way only a dispersion bistable response in the 
metamaterials is realized. A straightforward homogenization procedure is applied 
to describe such nonlinear systems as a composite medium with expression its 
material properties via some effective parameters. It should be noted that the main 
drawbacks of such structures are the low quality factor of resonances and considerable 
technological difficulties of their manufacturing in the optical range. 

Hence the metamaterials which are capable of supporting the trapped-mode 
resonant excitation [121 HOI HI] can be a successful alternative. The trapped-mode 
resonances appear in planar metamaterials designed on the basis of multi-element 
periodic arrays which typically consist of identical subwavelength metallic inclusions 
structured in the form of double-rings (DR's) j22j or asymmetrically split rings or 
squares [23l [25 • These particles are arranged periodically and placed on a thin 
dielectric substrate. The trapped-mode resonance is a result of the antiphase current 
oscillations in the particles of a periodic cell and they have high quality factor due 
to the weak interaction with free space. The nonlinearity can be included in such a 
system in the form of a substrate made of some Kerr-type or saturable absorption 
nonlinear material. 

This conception was recently confirmed experimentally in |25[ 126] . where the 
metamaterials which bear trapped-mode resonances were proposed to realize an 
enhancement of QD luminescence and electro-optic switching. The main advantage 
of the studied structures is the fact that a significant field enhancement is achieved 
in a thin planar system. Note that in these structures an active medium was used 
as a substrate of planar metamaterial. Also in our previous publications [27l [28] 
the features of a dispersion bistable response of metamaterials with asymmetrical 
and symmetrical particles placed on a nonlinear substrate were studied. It is shown 
that the effect of nonlinearity appears as the formation of closed loops of bistable 
transmission within the frequency of the trapped-mode resonance due to the strong 
field localization in the system. 

The goal of this paper is to show an enhancement of absorption bistability in 
a planar metamaterial which bears the trapped-mode resonances in sight to design 
the improved saturable absorbtion mirror. An array of metallic double-ring elements 
is chosen to construct the metamaterial due to the polarization insensitivity of this 
structure at the normal incidence of an exciting wave. 
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2. Problem statement and method of solution 

The square unit cell of the structure under study has a size d — dx — dy — 800 nm 
and consists of one metallic DR (Fig. [1]). The radii of the outer and inner rings are 
fixed at ai = 290 nm and 02 = 230 nm, respectively. The width of the both rings is 
2w = 40 nm. The array is placed on a dielectric substrate with permittivity e and 
thickness h = 150 nm. Suppose that the normally incident plane monochromatic wave 
has a frequency uj and a magnitude A. 



When the intensity of input light is small, i.e. in the linear regime, in order to 
calculate the optical response of metamaterial, two different techniques based on the 
method of moments (MoM) [29] and the pseudo-spectral time-domain (PSTD) method 
[23] were proposed earlier. The MoM is a numerical-analytical method that involves 
solving an integral equation which is related to the surface currents induced on the 
metallic pattern by the incident electromagnetic wave, then calculating the scattered 
fields produced by the currents as a superposition of partial spatial waves. In the 
framework of this method the metallic pattern is treated as a thin perfect conductor. 
On the other hand, the PSTD is a direct numerical scheme which uses a spatial 
spectral expansion of fields and the differentiation theorem for Fourier transforms to 
calculate the spatial derivatives. It is significant that the PSTD allows us to take 
into account the strong dissipation and dispersion of the metal permittivity of the 
elements. Nevertheless, in [24] it was revealed that the results obtained with the both 
techniques are in good agreement with each other down to the mid-IR region where 
the trapped-mode resonances are still well observed and have high Q factor. 

So we restrict our study within this region and use the MoM because it allows us 
to calculate the magnitude of currents J along the DR element since it is important 
to estimate the inner field magnitude inside the system. On the other hand, such 
a numerical-analytical treatment enables us to clearly understand the physics of the 
phenomena under investigation. Thus, in virtue of the MoM, the magnitude of the 
current J, the reflection i?, transmission T and absorbtion W = 1 — — jTp 
coefficients can be determined numerically as functions of frequency and permittivity 




Figure 1. (Color online) Fragment of the planar metamaterial of DR's and its 
elementary unit cell. 
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of substrate and presented in the symbolic form 

u^Fu{u:,e), (1) 

where u — J, R, T, or W. 

Next we should note that, at the trapped-mode resonance, due to specific current 
distribution on the metallic pattern, the electromagnetic energy is symmetrically 
confined to a very small region between two half-arcs of the rings, where the energy 
density reaches substantially high values. The latter feature allows us to construct 
approximate model based on the transmission line theory to estimate the inner field 
intensity caused by rings 28 . According to this theory, in each unit cell, the DR 
is symmetrically divided into two parts where the half-arc of the inner and outer 
conductive rings are considered as two wires with a distance 6 = ai — 02 — 2w between 
them. Along these wires the currents with equal magnitude flow in opposite directions. 
Thus the electric field strength is defined as 

E^n = V/b = ZJ/b, (2) 

where V is the line voltage, J is the magnitude of current induced in the DR element, 
and Z is the impedance of line. The impedance is determined at the resonant frequency 
aeo — d/Xo as, 

_ dOleep 
dCo ' 

where / — 7r(ai + a2)/2, and 



Co^iln 



2w V \2wJ 



is the capacity in free space per unit length of line, p = ai — 02 . In view of the small size 
of the translation cell of the array in comparison with the wavelength, we suppose that 
the nonlinear substrate remains to be a homogeneous dielectric slab under intensive 
light. Thus, from this model it is possible to evaluate the relation between the inner 
electric field strength and the current magnitude J averaged along one half-arc of the 
DR-element. 

Next, in the nonlinear regime^ we suppose that the structure's substrate is a 
material whose absorption property depends on the intensity of the electric field 
lin = (l/2)7T.eoc|£^iTiP inside it. Here the complex refractive index n = -Je of the 
substrate is defined in the form [3^ 

z c 

n = n' ^ in" = n' + -z—ct, (3) 
2 u) 

where a is an absorption index. Thus the permittivity e of the substrate can be found 
as a function of a, e = e(a). 

As a saturable absorbing material for the structure's substrate we consider 
chalcogenide glasses doped with quantum dots of other elements such as Ge, As and 
Sb due to their wide availability for the mid-IR range (HI E]- It is known that the 
systems based on quantum dots have discrete energy level. To describe the effects of 
saturation in such systems a two-level model is widely used. Further, the simple two- 
level model is developed in three directions. The first direction consists in use of the 
three- and four-level models, which is important in studying the process of creating 
a population inversion in lasers. The second direction of the model development is 
taking into account the inhomogeneous broadening of spectral lines due to the fact that 



Enhancement of absorption bistability by trapping light planar metamaterial 



6 



every quantum structure element (molecule, ion, quantum dot, etc.) is in the fields 
of its neighbors, which may be different. Finally the third direction is considering 
the dynamic Stark shift of the energy levels of molecules under the action of light. 
As usual, the two-level model is preferable when studying bistable devices because in 
many cases it gives good agreement with experiments [6l [7] and due to its simplicity. 

Thus to describe the dependence of permittivity e on the absorption index a we 
use a two-level model of a fast relaxing absorber whose absorption index a is described 
by the formulae |31] 



0.25ao7^ 

- Wo)2 +0.257^(1 + UnlLat) 

where ao is the weak-field absorption index, 7 is the weak-field absorption line 
broadening, l^at is the saturation intensity and is the frequency of the center 
of the spectral absorption line. The absorption curve defined by equation ([5]) has a 
Lorentzian shape, and when /i„ = l^^^t the absorption index is reduced to half of the 
peak of the absorption line. Here the terms describing the dynamic Stark shift are 
discarded under an assumption that the effect of the dynamic Stark shift occurs under 
light intensities, much larger than those required for saturation of absorption ^31 . The 
absorption line broadening depends on the inner field intensity as 7-y/l -I- 1^/ Isat- In 
the center of the spectral line, the absorption index is a{ijJo,Iin) = oto/(\ + Im/Isat), 
whereas on the distant sides of the absorption line, oj — wo 3> 7, the value of a{Lo, Im) 
is practically the same as a(w,0), i.e., under a certain intensity level, the most 
pronounced decreasing of absorption is reached in the center of the line. Hence 
it follows that the best way to obtain strong interaction of the metamaterial with 
a nonlinear substrate lies in tuning the trapped-mode resonant frequency to the 
frequency of the center of the spectral line of a saturable material. This fact is 
also proved with experimental results of the QD's luminescence enhancement in the 
plasmonic metamaterial |25| . 

As the absorption index a is a function of the inner field intensity /i„, so, too, 
is the permittivity of substrate, e = e{Ii„). According to equation ([2|), the inner field 
intensity is a function of the averaged current magnitude, lin = Iin{J)- Thus the 
nonlinear equation on the averaged current magnitude in the metallic pattern can be 
obtained in the form [ 27, 

J = i-Fj(w,£(/,„(J))), (5) 

where A is a dimensionless coefficient which depicts how many times the incident 
field magnitude A is greater than 1 V cm~^. Thus, the magnitude A is a parameter of 
equation ([5]), and, at a fixed frequency w, the solution of this equation is the average 
current magnitude J which depends on the magnitude of the incident field (J = J{A)). 

Further, on the basis of the current J (A) found by a numerical solution of equation 
(O, the permittivity of the nonlinear substrate e = e{Iin{A)) is obtained and the 
reflection, transmission and absorption coefficients are calculated as functions of the 
frequency and magnitude of the incident field. 



3. Numerical results and discussion 



Typical frequency dependences of the current magnitude, reflection, transmission 
and absorption coefficients calculated in the linear regime under the assumption of 
a nondispersive substrate with the method of moments are given in figure. [2] One 
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Figure 2. Frequency dependences (ae = d/X) of magnitude squares of the average 
current (in a.u.) (a), transmission (b), reflection (c) coefficients, and absorption 
coefficient (d) in the case of linear regime at a low intensity of the incident field. 



can see that a sharp resonance appears in the band where the maximum of current 
magnitude occurs. At this frequency band, the induced currents in the inner and 
outer rings oscillate in opposite directions. An electromagnetic field of two closed 
strips of the rings is similar to the nonradiating field of a double- wire line of resonant 
length. Thus, the oppositely directed but almost equal currents of the inner and 
outer rings of the array yield an electromagnetically trapped mode. The scattered 
field produced by such a configuration of current is very weak in comparison with 
an intensity of the stored field, and, as a consequence, the coupling of the planar 
metamaterial to free space is small and, therefore, radiation losses are reduced, which 
ensures a high quality factor resonant response. Evidently, as the value of Joule losses 
in the metamaterial increases, the magnitude of current and quality factor of the 
resonance decrease, but, nevertheless, the resonance remains to be well observed |32j . 
By this means, we suppose that, near the frequency of trapped-mode excitation, such 
a metamaterial will strongly react on the intense incident wave when the substrate 
is made of a nonlinear dielectric even in the case when the dielectric is dissipative. 
Also, in figure[2rb) (figure[Hc)), one can see the evolution of approximately symmetric 
close to Lorentzian profile of the transmission (reflection) resonance to an asymmetric 
Fano one [Ml EH ES] when the losses decrease which is very suitable to realize abrupt 
switching in the nonlinear regime. 

Thus, in the case of the array placed on the substrate whose complex permittivity 
depends on the field intensity, the dependences of the inner intensity versus the 
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Figure 3. (Color online) Dependences of the inner intensity versus the incident 
field magnitude in the case of the nonlinear permittivity of the substrate. For this 
and further calculations, we use the experimental results obtained in [B]. Thus, 
the saturation intensity is chosen to be Isat = 60 kW/cm^, which is typical for the 
phosphate glasses doped with 6-nm PbS QD's in the IR-range, 7 = 7(i/(27rc) = 1, 
n' = \/4.1. The dimensionless frequency ae corresponds to the trapped-mode 
resonance and to the peak of the absorption line, se = aeo = uiod/ (2itc) = 0.3. 



incident field magnitude /.„ = /m(A) have a typical form of hysteresis loops (figure [3]). 
The origin of the absorption bistability is well studied in the theory of Fabry-Perot 
resonators [1], which is also applicable to the nonlinear metamaterial. In accordance 
with this theory, the nature of the hysteresis curve can be cleared from the next 
considerations. A slight increase in the incident field intensity causes a slight decrease 
in the metamaterial substrate absorption which, in turn, permits the pattern to accept 
more power. This added power further saturates the saturable absorber which again 
permits more power to be coupled in the pattern. This process is cumulative, with 
the result that the change in the power level in the system exceeds the original small 
change in the incident intensity. Thus, there is a threshold level at which the system 
is unstable and abruptly switches states. A similar process is involved as the input 
intensity is reduced, causing the system to switch abruptly from a low attenuation 
state to a high attenuation state at some lower threshold level. 

From [4] the condition of achieving completely absorption bistability is also 
known, a^L/T > 8, where L and T are the length and transmission coefficient of the 
Fabry-Perot resonator, respectively. From this condition it follows that the switching 
threshold essentially depends on the value of the absorption index ag. Thus, from 
our calculation (figure [3]) it is clear that, for the chosen metamaterial parameters, the 
switching appears only when ao is sufficiently high since the substrate is thin. 

Nevertheless there is an important distinctive feature of the studied structure in 
contrast to the Fabry-Perot one. This peculiarity is that the resonance in the Fabry- 
Perot cavity has the Lorenzian-shape whereas the trapped-mode resonance in the DR- 
structure has the Fano-shape, and the latter one is characterized by an asymmetrical 
peak-and-trough profile. 

From the classical point of view [33], the resonant conditions of the studied 
structure can be described using two weakly coupled harmonic oscillators, where one 
of them is driven by a periodic force (see, for an example, figure 2(a) in [34 ). In such 
a system, in general, the spectrum of the forced oscillator consists of two resonances. 



Enhancement of absorption bistability by trapping light planar metamaterial 



9 



and while the first resonance is characterized by a symmetric profile, described by a 
Lorenzian function, the second one is characterized by an asymmetrical Fano profile. 
At a certain frequency, the amplitude of the forced oscillator becomes zero, as a result 
of destructive interference of oscillations from the driving force and coupled oscillator. 
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Figure 4. Dependences of the magnitudes of the reflection, transmission and 
absorption coefficients versus the incident field amplitude for different absorption 
index oq in tliG ca,SG of the nonliiiGcir permittivity of the siibstra-tci Isat — 
60 kW/cm^, 7 = 7d/(27rc) = 1, n' = V4T, eg = aeo = 0.3. 

In our case, the different current distributions in the metallic pattern lead to the 
appearance of similar resonant conditions (figure [2jb)). Thus, the resonant feature I is 
a result of resonant excitation of a predominantly dipole mode (i.e., current oscillations 
symmetric with respect to the symmetry axis of DR) in either outer or inner ring 
and this resonance is described by a Lorenzian function. At the feature II, both 
rings are excited equally, while the induced currents in the inner and outer rings 
oscillate in opposite phases, yielding an electromagnetically trapped mode [22j |32] 
whose resonance has an asymmetrical profile. 

Such asymmetrical spectral profile of the transmission coefficient which vary from 
low to high over a very narrow frequency range can be very useful in all-optical 
switchings since there are gently sloping bands of the high reflection and transmission 
before and after the resonant frequency [34j [35] . This peculiarity of spectra manifests 
itself in the stepwise variation of magnitudes of reflection, transmission and absorption 
coefficients calculated versus the incident field amplitude (figure H} . These figures 
also illustrate another effect of the shape asymmetry of the trapped-mode resonance, 
namely the formation of the hysteresis with closed loops in the curves of both reflection 
and absorption coefficients (figure |4}d). In particular, closed loop appears when the 
level of refiection (absorption) before the threshold amplitude is greater then that one 
after it. Thus when the hysteresis is formed, the part of the curve that corresponds 
to the rising edge remains above the part of the curve that corresponds to the falling 
edge. 

4. Conclusions 

In conclusion, a planar nonlinear DR-metamaterial, which bears a high quality factor 
Fano-shape trapped-mode resonance is a promising object for a realization of a 
polarization insensitive absorption bistable switching. The main advantage of the 



Enhancement of absorption bistability by trapping light planar metamaterial 



10 



studied structure is the possibility to achieve the bistable transmission at low input 
intensity, due to a large quality factor of the trapped-mode resonance. 
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